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ABSTRACT 

Advancements in infrared (IR) interferometry open up the possibility to spatially resolve active 
galactic nuclei (AGN) on the parsec-scale level and study the circumnuclear dust distribution, com- 
monly referred to as the "dust torus", that is held responsible for the type 1/type 2 dichotomy of 
AGN. We used the mid-IR beam combiner MIDI together with the 8 m telescopes at the Very Large 
Telescope Interferometer (VLTI) to observe the nucleus of the Seyfert 2 galaxy NGC 424, achieving 
an almost complete coverage of the uw-plane accessible by the available telescope configurations. We 
detect extended mid-IR emission with a relatively baseline- and model-independent mid-IR half-light 
radius of (2.0 ± 0.2) pc x (1.5 ± 0.3) pc (averaged over the 8 — 13 ^m wavelength range). The extended 
mid-IR source shows an increasing size with wavelength. These properties are in agreement with the 
idea of dust heated in thermal equilibrium with the AGN. The orientation of the major axis in position 
angle ~ —27° is closely aligned with the system axis as set by optical polarization observations. Torus 
models typically favor extension along the mid-plane at mid-IR wavelengths instead. Therefore, we 
conclude that the majority of the pc-scale mid-IR emission (>60%) in this type 2 AGN originates 
from optically-thin dust in the polar region of the AGN, a scenario consistent with the near- to far-IR 
SED. We suggest that a radiatively-driven dusty wind, possibly launched in a puffed-up region of the 
inner hot part of the torus, is responsible for the polar dust. In this picture, the torus dominates the 
near-IR emission up to about 5 /xm, while the polar dust is the main contributor to the mid-IR flux. 
Our results of NGC 424 are consistent with recent observations of the AGN in the Circinus galaxy 
and resemble large-scale characteristics of other objects. If our results reflect a general property of 
the AGN population, the current paradigm for interpreting and modeling the IR emission of AGN 
have to be revised. 

Subject headings: galaxies: active - galaxies: Seyfert - galaxies: individual: NGC424 - infrared: 
galaxies - techniques: high angular resolution 



1. INTRODUCTION 

The dusty environment around supermassive black 
holes in active galactic nuclei (AGN) came in reach 
of direct observations using the new high spatial res- 
olution capabilities of long-baseline infrared (IR) in- 
terferometry. Within the last years several observing 
campaigns revealed emission sources in the nuclear re- 
gion that shows clear characteristics of AGN-heated 
dust on the sub -parsec to parsec scale in the near- 
and mid-IR (e.g. Jaffc et al. '2004; 'Tristram et al.' '2007; 
Beckert et al. 2008; Kishini oto ct aL .2009a; Raba n et al.„ 
20091: iTristram et' all [20091: TBurtscher et al.|[2009[) . The 
two prototypical type 2 AGN in NGC 1068 and the 
Circinus galaxy have been studied in detail using the 
mid-IR beam combiner MIDI at the Very Large Tele- 
scope Interferometer (VLTI) and pairwise combinations 
of the four 8 m-telescopes. On these rather long baselines 
both objects were highly resolved with rather low visibil- 
ities that have been interpreted as the reemission signa- 
tures of the circumnuclear obscuring region, commonly 
dubbed "dust torus". This geometrically-thick torus is 
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an essential part of the unification scheme of AGN and 
explains the difference of type 1 and type 2 AGN by 
angle-dependent obscuration. Recently, IR interferome- 
try has shown that the spectral and spatial characteris- 
tics of the mid-IR emission sources (spectral slope and 
apparent size) in type 1 AGN are tightly correlated with 
the brightness profile and radial distribution of the dust 
(jKishimoto ct al. 2011b). 

One problem with these previous observations is that 
they either miss most of the single-aperture flux, which 
has been resolved out (NGC 1068 and Circinus), or 
that they do not have the position-angle coverage nec- 
essary to constrain the shape/orientation of the nu- 
clear emission source. This structural information is 
an important input for torus models that simulate IR 
images and predict certain position-angle dependent 
spectral and spatial feat ures (jSchartmann et al] 120081 : 
IHonig fc Kishimotol[20T0l ). 

In this paper, we present mid-IR interferometry ob- 
servations of the Seyfert 2 galaxy NGC 424. These 
observations cover almost the entire uu-plane accessi- 
ble with the VLTI UTs. NGC 424 is a SBO/a galaxy 
at a distance of 45.7 Mpc. The galactic disk is in- 
c lined at about 70° bas ed on an axis ratio of 0.36 
(Kirhakos fc S teineii[T990f ) . It hosts an AGN that is clas- 
sified as an optical type 2 A GN with polarized broad lines 
(i.e., a hidden type 1 AGN: iVeron-Cettv fc Veronll2010D 
based on spectro-polarimetric observations (jMoran et al.l 
120001 ). The gas column towar ds the nucleus is quite 
large: ICoUinge fc Brandtl (|2000D report that the Hydro- 
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Figure 1. iii;-plane coverage (on-sky orientation) of the mid-IR 
intcrfcrometric observations of NGC 424. North is up, east is left. 



gen column density is probably in the Compton- 
thick regime (> 2 X 10^4 cm-2), while E amassa et al.l 
()2011|) find less, but still significant. X-ray obscuration 
[Nh ~ 1.7 X 10^'^ cm~^). Both values are consistent with 
a heavily obscured AGN, and most of the obscuring col- 
umn is probably intrinsic to the AGN, in spite of the con- 
siderable inclination of the h ost galaxy. Based o n the ve- 
locity dispersion of the bulge, iBian fc Gul (|2007f l estimate 
the mass of the central black hole as logMen = 7.78, 
leading to an Eddington ratio of /Edd = 0.13, which is 
typical for a Seyfert gal axy. The nucleus is radio-quiet 
and unresolved in radio (jMundell et al.|[200ClD . 

In Sect. [21 we present the interferometric data and de- 
scribe our data reduction strategy. For interpretation we 
also compiled a high spatial resolution IR SED of single- 
telescope data to complement the interferometry. In 
Sect. [21 wavelength- and position- angle-dependent sizes 
are extracted from the interferometric observations and 
the IR SED is analyzed. Both interferometry and pho- 
tometry data are simultaneously modeled in Sect.[4l The 
results are discussed in Sect. [5] and put into context of 
AGN unification. Finally, we summarize our findings in 
Sect, [i 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Interferometry 

The interferometric observations have been carried out 
during three different campaigns in 2007, 2009, and 2010 
using the mid-IR beam combiner MIDI at the Very Large 
Telescope Interferometer (VLTI) in Chile. This instru- 
ment allows for spectro-interferometric observations in 
the mid-IR waveband from 8 — 13 /im. We used sev- 
eral different combinations of the four 8-m telescopes to 
achieve various position angles and projected baseline 
points to fill the accessible wti-plane of NGC 424. In 



Fig. [T] we show the ww-plane coverage of our observa- 
tions. The uu-plane is shown as projected on the sky to 
allow for an easy comparison with position angles. While 
the shorter baselines of up to 60 m cover all quadrants, 
the longer baselines are oriented preferentially in NE/SW 
direction, in line with the location of the UT telescopes 
at the observatory. As illustrated, our observations can 
be considered complete for this target considering the 
achievable coverage using the VLTI UTs. 

The data have been extracted by the standard 
MIA+EWS package developed for MIDI. Reduction and 
calibration were perforn ied using the SNR maxim ization 
method as described in iKishimoto et al.l ()2011bf) . Each 
individual fringe-track data set and a dedicated inter- 
ferometric calibrator star have been reduced together to 
determine the required smoothing of the calibrator data 
to achieve the best signal-to-noise ratio of the correlated 
fluxes. During all three observing campaigns, we ob- 
served the same calibrators near NGC 424 (HD9362 and 
HD8498) for consistency, except for one set in 2007 where 
the only available calibrator was observed at significant 
spatial and temporal distance from the science target 
(HD219449). Cross-calibration showed that HD9362 and 
IID8498 remained essentially constant in flux and size 
over the observed time span. 

When two or three science fringe track sets have been 
observed together (i.e. within a 15min interval), they 
were averaged after data reduction with weights accord- 
ing to their individual SNR. This led to 3 independent 
correlated flux measurements in 2007, 17 in 2009, and 
27 in 2010 as listed in Table [H The SNR values of 
all data sets range from 3.0 to 7.2 with a median of 
4.3. This is slightly hig her than for a typical AGN (see 
IKishimoto et al.l2011bf ) and can be explained by the mid- 
IR brightness of the source. 

We extracted the total fluxes from all of the photom- 
etry data that has been taken together with the fringe 
tracks. As sky reference we used the background flux 
from a 5 pixel offset position to the peak of the tar- 
get flux. Calibration was done using the same calibrator 
star as used for the interferometric data. All data sets 
of each year have been averaged to obtain a reference 
total flux of NGC 424 in the 8-13 ^J.m range for that 
year. The 2007, 2009, and 2010 spectra are shown in 
Fig. [2] While the 2009 and 2010 fluxes are consistent, 
the 2007 spectrum is systematically lower by a factor of 
0.8 mostly independent of wavelength. 20-30% variabil- 
ity over a span of 2 years may be possible, but given 
the bin-to-bin (systematic) scatter in the 2007 observa- 
tions, it is more likely an observational inaccuracy. In- 
deed, the 2007 fluxes are based on only 3 photometry 
data sets while the 2009 and 2010 spectra are averages 
over 23 and 22 sets, respectively. For a good total flux 
reference, we calculate a "master" total flux spectrum, 
i^tot, from the 2009 and 2010 data by averaging both 
sets and spline-fitting the atmospheric ozone feature at 
around 9.6 /im. This spectrum will be used as the ref- 
erence for all correlated fluxes Fcorr;i- The visibilities 
that we will use in the data analysis are then computed 
as Vi — Fcori-i/ Ftot, and visibility errors are propagated 
based on the total and correlated flux uncertainties. As 
described in IKishimoto et al.l ()2011bD , these uncertain- 
ties contain the fluctuations of the individual correlated 
and total flux frames, the differences to the mean when 
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Table 1 

Observation log and data properties 



Observing Date 


telescopes 


PBL 


PA 


calib. 


max. 




(UT) 








(m) 


(deg) 




SNR 


2007- 


-10- 


-24 01:53 


U3- 


-U4 


56.1 


94.2 


HD9362 


4.1 


2007- 


-10- 


-24 03:57 


U3- 


-U4 


63.0 


112.3 


HD9362 


6.2 


2007- 


-10- 


-25 04:00 


U2- 


-U3 


45.5 


43.3 


HD219449 


3.5 


2009- 


-08- 


-02 


06:37 


U3- 


-U4 


57.7 


88.0 


HD8498 


4.4 


2009- 


-08- 


-02 


07:37 


U3- 


-U4 


57.8 


97.1 


HD8498 


5.1 


2009- 


-08- 


-02 


08:48 


U3- 


-U4 


61.8 


107.4 


HD9362 


6.0 


2009- 


-08- 


-02 


09:42 


U3- 


-U4 


62.4 


115.5 


HD9362 


5.0 


2009- 


-08- 


-03 


06:30 


Ul- 


-U2 


56.5 


6.3 


HD8498 


5.4 


2009- 


-08- 


-03 


07:19 


Ul- 


-U2 


56.3 


13.1 


HD8498 


6.5 


2009- 


-08- 


-03 


08:26 


Ul- 


-U2 


55.9 


22.0 


HD9362 


6.8 


2009- 


-08- 


-03 


09:50 


Ul- 


-U2 


54.2 


31.7 


HD9362 


7.2 


2009- 


-08- 


-04 


06:23 


Ul- 


-U3 


102.4 


11.0 


HD8498 


3.0 


2009- 


-08- 


-04 


07:37 


Ul- 


-U3 


102.2 


22.0 


HD9362 


3.0 


2009- 


-08- 


-04 


08:52 


Ul- 


-U3 


100.7 


32.3 


HD9362 


3.2 


2009- 


-08- 


-04 


09:46 


Ul- 


-U3 


98.2 


38.3 


HD9362 


3.4 


2009- 


-08- 


-05 


05:39 


Ul- 


-U4 


122.0 


25.7 


HD8498 


4.0 


2009- 


-08- 


-05 


06:31 


Ul- 


-U4 


125.7 


36.5 


HD8498 


3.8 


2009- 


-08- 


-05 


07:23 


Ul- 


-U4 


128.8 


46.2 


HD8498 


3.9 


2009- 


-08- 


-05 


08:30 


Ul- 


-U4 


130.1 


57.2 


HD9362 


4.1 


2009- 


-08- 


-05 


09:08 


Ul- 


-U4 


128.7 


62.7 


HD9362 


4.2 


2010- 


-10- 


-19 


01:31 


Ul- 


-U3 


102.4 


11.9 


HD9362 


3.6 


2010- 


-10- 


-19 


02:35 


Ul- 


-U3 


102.3 


21.4 


HD9362 


3.5 


2010- 


-10- 


-19 


06:11 


U3- 


-U4 


59.8 


131.9 


HD8498 


4.4 


2010- 


-10- 


-19 


07:08 


U3- 


-U4 


56.9 


143.9 


HD8498 


3.9 


2010- 


-10- 


-19 


08:02 


U3- 


-U4 


54.2 


157.1 


HD8498 


3.8 


2010- 


-10- 


-20 


01:06 


U2- 


-U4 


78.1 


50.5 


HD8498 


4.8 


2010- 


-10- 


-20 


01:58 


U2- 


-U4 


83.7 


61.3 


HD9362 


3.8 


2010- 


-10- 


-20 


02:50 


U2- 


-U4 


87.7 


70.5 


HD9362 


4.4 


2010- 


-10- 


-20 


03:49 


U2- 


-U4 


89.4 


79.9 


HD8498 


4.3 


2010- 


-10- 


-20 


04:40 


U2- 


-U4 


88.0 


87.8 


HD9362 


4.3 


2010- 


-10- 


-20 


05:32 


U2- 


-U4 


83.9 


95.9 


HD8498 


4.5 


2010- 


-10- 


-20 


06:17 


U2- 


-U4 


78.2 


103.7 


HD9362 


3.8 


2010- 


-11- 


-28 


00:47 


Ul- 


-U4 


130.2 


55.4 


HD8498 


4.7 


2010- 


-11- 


-28 


01:38 


Ul- 


-U4 


128.7 


62.9 


HD8498 


4.8 


2010- 


-11- 


-28 


02:30 


Ul- 


-U4 


123.6 


70.0 


HD8498 


4.8 


2010- 


-11- 


-28 


03:00 


Ul- 


-U4 


119.0 


73.9 


HD9362 


4.6 


2010- 


-11- 


-28 


03:39 


Ul- 


-U4 


111.2 


79.0 


HD9362 


4.4 



Table 2 

Nuclear IR photometry of NGC 424 extracted from archival 
VLT/ISAAC images. 



filter 


A 


flux 


calibrator 


obs. date 




(^tm) 


(Jy) 






(UT) 


Js 


1.24 


0.0043±0.0003 


S754-C 


2002- 


-11-14 01:01 


H 


1.65 


0.0177±0.0007 


S754-C 


2002- 


-11-14 01:05 


Ks 


2.16 


0.0445±0.0011 


S754-C 


2002- 


-11-14 01:10 


L 


3.78 


0.2132±0.0046 


HD2811 


2003- 


-08-17 07:12 


M 


4.66 


0.3038±0.0294 


HD2811 


2003- 


-08-17 07:48 



Note. — The fractional contributions of the nuclear flux to the 
total flux in the central I'.'O are 0.36, (Js-band), 0.57 (H-band), 
and 0.65 (Ks-hand) based on the nuclear photometry after 2-D 
host subtraction with a Scrsic profile in each band. 



SNR-weighting several data sets, and a 5% systematic 
component resulting from the analysis of low-SNR data 
(see Appendix Fig. A. 5 in iKishimoto et 3111201 Ibl ). 

2.2. Photometry 

In addition to our MIDI interferometry, we collected 
photometric data mainly from the ESO archive and lit- 
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Figure 2. Total 8 — 13 fira flux spectrum of NGC 424 as measured 
with MIDI. The red-dashed line represents the data from 2010, 
the blue-dashed line from 2009, and the green-dash-dotted line 
shows the 2007 data. The black-solid line is the error-weighted 
average of the 2009 and 2010 data and used as the reference for 
the interferometric data. The spectral feature around 9.6 fim in 
the 2009 and 2019 data is a result of non-perfect subtraction of the 
atmospheric ozone feature. 

erature to build a generic IR SED of the nuclear emis- 
sion of NGC 424. The ESO data archive holds near-IR 
imaging data of NGC 424 in the Js-, H-, Ks-, L- and 
Af-bands observed with ISAAC. We reduced the data us- 
ing standard procedures. For the three near-IR bands, 
we first subtracted the host galaxy from the nucleus by 
fitting a Sersic profile to the extended (> 1") component 
(i.e. masking the nuclear region for the fit) before ex- 
tracting nuclear fiuxes. The resulting fluxes are shown 
in Table. H 

It is important to point out that the L- and Af-band 
data from the 8m telescope with a (seeing-limited) spatial 
resolution of about 0'.'4 are consistent, within error, with 
the 3.6 /ini and 4.5 ^m Spitzer data with a spatial reso - 
lution of about 1'.'2 reported bv lGallimore et all ()2010D . 
This can be taken as evidence that the IR emission in 
the central arcsecond is dominated by emission from the 
AGN without significant enhancement by starformation 
in the nuclear region. Indeed, the Spitzer IRS spectrum 
of the nuclear ^ 4" shown in iGallimore et al.l (|201Q) dis- 
plays PAH features with only small equivalent widths, 
the large st being the 7.7 jum feature but no line visible at 
ll.S/im ()Wu et al.ll20T0l ; AGN fraction > 0.9 using the 
authors' relation between AGN contribution, IRAC and 
IRAS photometry). 

3. RESULTS 

3.1. Broad-band infrared SED 

In Fig. [3] we show a broad-band IR SED from about 
1 i xm. to 90 fim combinin g the Spitzer data as shown 
in IGallimore et al.l ()2010D with the 8m-telescope data 
presented in this work. The plot illustrates the IR 
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Figure 3. 1 - 90 broad- band SED of NGC 424. The gray lines 
(IRS) and ope n circles (IRAC a nd MIPS) show the Spitzer data 
as presented bv lGallimore et al.l (^OlO) with a spatial resolution of 
> 1.2" on the shorter and > 10" on the longer wavelength side. 
The red-filled circles represent the VLT/ISAAC data with a spatial 
resolution of ~ 0"4 while the blue line is our VLTI/MIDI total fiux 
reference extracted from a window of about 0"5 X (y.'7. The years 
in which the observations were taken are indicated. 



bump typical of dust emission and seen in all radio- 
quiet AGN. According to the unification scheme, it is 
caused by circumnuclear dust in the equatorial region 
of the AGN that re-emits the absorbed "big blue bump" 
(BBB) UV/optical radiation (presumed to originate from 
the accretion disk) . One of the interesting properties of 
this target is that we see the IR SED rising at short 
"wavelengths, flattening in the 8 — 20 /im region, and then 
falling off towards 100 ^m, despite the comparable large 
intrinsic apertures involved. The resolution of the MIPS 
data at 70 /urn corresponds to '^4.6 kpc. Yet the SED 
is more or less a downward continuation of the data at 
shorter wavelengths and much higher angular resolution 
(at least to first order; see below). Indeed, the whole 
SED from the near- to far-IR connects smoothly despite 
the different instruments and resolutions involved. We 
interpret this as evidence that the unresolved IR nucleus 
dominates all of the IR emission without any significant 
contribution from other sources, e.g. circumnuclear star- 
forming regionfQ, and that the AGN is the only source 
of energy that produces the IR emission in the nuclear 
region. 

There are two notable exceptions from the smooth con- 
nection of all SED data. The first one is a slight offset of 
the MIPS data with respect to the longest wavelengths 
in the IRS data. This may be an indication that there 
exists an additional, cool IR component contributing to 
longer wavelengths. If this were true, then we would ex- 
pect that the slope changes from the downward tail of 

^ If there were significant starformation within the photometric 
aperture, we would expect both strong PAH emission features and 
a bump in the far-IR. Neither is seen here though. 



the IRS data through the MIPS data. However, this is 
not the case. In fact the data seem to be just offset by a 
factor of 1.3, which could be a sign of an issue with flux 
calibration in the far-IR, but we cannot exclude the pos- 
sibility of slightly extended emission between the smaller 
and larger apertures. Anyway, the flux is strongly and 
monotonically decreasing toward the far-IR. 

The second discontinuity in the SED is more interest- 
ing. While the ISAAC L- and Af-band measurements are 
well in agreement with the IRAC data and smoothly con- 
nect to the IRS short wavelengths, the MIDI total flux 
spectrum is significantly lower than the IRS data. The 
difference of a factor of 1.3-1.4 is difficult to be attributed 
to calibration problems, since both the IRS and MIDI 
calibration seems to be consistent in itself (for MIDI see 
Sec. 12. 1|) . Under the assumption that there are no un- 
known systematic problems with the MIDI data, we can 
imagine two possibilities for this discrepancy: first, there 
may be extended warm emission outside of the MIDI 
aperture of ^ 0'.'6 or 130 pc. The dust cannot be too hot 
(>800K) or too cold (<200K), however, because we do 
not see a disagreement between small and large apertures 
in the near-IR and also no bump in the IRS spectrum be- 
yond 20 /im. Moreover, the SED in the difference spec- 
trum is fiat, so that we can estimate a temperature of the 
possible extended emissio n of^ 300 — 400 K. When using 
the rsub — ^12 /im-relation (i Honig et al.ll2010l) to estimate 
the dust sublimation radius for this type 2 source, we ob- 
tain about 0.065 pc (with significant error bars). Given 
the constraints from the MIDI and Spitzer aperture, the 
emission responsible for the mid-IR fiux difference would 
then be located at a distance > 2000 rgub- When assum- 
ing a standard ISM dust composition, we would expect 
an equilibrium temperature with the AGN radiation of 
<100K at that distance, which contradicts the require- 
ments based on the SED shape as just discussed. Finally, 
an additional heating source that provides energy for the 
dust in a very narrow wavelength range is also difficult 
to imagine, in particular when required not to leave any 
other traces. We repeat that starformation is a very un- 
likely candidate because of the almost complete lack of 
PAH features (except a small feature at 7.7 /im) and the 
lack of a far-IR bump. 

The second possibility is mid-IR variability. In Fig. [3] 
we indicated the years in which each of the sets of obser- 
vations were taken. Interestingly the MIDI data is sepa- 
rated by about 4 years from the other observations. We 
recently discussed how IR variability of a dusty structure 
comp ares at different wavelengths (|H6nig fc Kishimotol 
1201 If ). The principle of IR variability is simple: if the in- 
cident radiation onto the dust (here: the BBB radiation) 
changes, the dust temperature at a given distance from 
the source changes as well and, consequently, the fiux 
varies. The main difference between near- and mid-IR 
wavelengths is the size and distance of the emitting re- 
gions. The near-IR emission is usually confined to a very 
narrow region close to the sublimation radius. Therefore 
the light-crossing time of this region for any variability 
signal is short and the initial signal is not significantly 
diluted {— comparable amplitudes of initial BBB signal 
and near-IR signal). The mid-IR emission is more com- 
plex. It is a combination of the Rayleigh- Jeans tail of 
hotter dust and the peak in emission from cooler dust. 
Moreover, the radial distribution of the dust also mat- 
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Figure 4. Comparison of the 8 fira normalized mid-IR SED of 
NGC 424 (black solid line) to AGN templates. The red-dashed 
line is an AGN template from MuUanev et al. (2011) combining 
type 1 and type 2 AGN . The light-blue lines are 10 nearby AGN 
from IHonig et al.l |[20TOh observed with VLT/VISIR. These data 
have an intrinsic resolution between 30 pc and 80 pc at 8 ^m to 
approximately match the resolution of NGC 424 (55 pc). The dark- 
blue dash-dotted line is the median spectrum of these AGN. 



ters. Overall the mid-IR emission originates from a rela- 
tively large range of distances from the AGN, so that the 
variability signal is smeared out and the amplitudes are 
smaller, unless the initial signal has a very long duration 
(e.g. a long-ter m break or jump in the A GN luminos- 
ity; see Fig. 2 in IHonig fc Kishimotol[20lll for an exam- 
ple). Therefore, mid-IR variability is a viable possibility 
if both (1) the amplitude of variability and (2) the size of 
the region wher e this v ariability originates from are phys- 
ically plausible. iGlassI p004 ) reported near-IR variability 
in the K-hand: the host-uncorrected K-hand magnitude 
showed peak-to- valley variation of 0.34 mag over several 
years. Comparing the uncorrected K-hand magnitude 
with our nuclear i^-band flu x rev eals that about 70% 
of the iiT-band flux of Glasi (|2004[ ) originates from the 
host. The corrected K-hand variability is, therefore, a 
factor of ~ 1.8 using our nuclear if-band measurement. 
This near-IR amplitude is sufhcient to explain a factor 
of 1.4 variability in the mid-IR. The second necessary 
condition related to the size of the mid-IR emission re- 
gion can be constrained by the elapsed time between the 
epochs of the mid-IR observations: Within the 4 years, 
the light propagated about 1.2 pc, so that in order to 
obtain the potential decrease in mid-IR flux by a factor 
of 1.4, the bulk of the mid-IR emission must come from 
within this radius. We can test this possibility with our 
interferometric data (see below). 

The SED of NGC 424 is relatively blue when compared 
to other type 2 or even type 1 AGN. In Fig. Uwe show 
the mid-IR emission of NGC 424 (black solid line) and 
two different AGN templa te SEDs. The re d-dashed line 
is an AGN template from iMuUanev et al.l (|201 h ) where 



Spitzer IRS data of 25 AGN were corrected for galac- 
tic/starformation and averaged without distinguishing 
among types. We already showed that the Spitzer IRS 
flux of NGC 424 is consistent with the IR fluxes from 
ground-based 8 m data. These data have higher angu- 
lar resolution than Spitzer data. At 8 /im, the PSF in 
NGC 424 has a size of ^5 5 pc. To account for possi- 
ble resolution effects in the iMullanev et a l.l (12011') tem- 
plate, we used the 10 AGN from lHonigeTa l. (2010, light- 
blue lines) with similar resolution as NGC 424 (range 
of 30— SOpc at 8 /im) and calculated a median template 
(dark-blue dash-dotted line). NGC 424 shows a bluer 
SED than the combined SED. The spectral index a (de- 
fined as oc i^") is a ^ —2.1 for type 2s a nd a ~ 
for ty pe Is, approximately in agr eement in IHonig et al.l 
(pom ) and MullanejLeLaU (poTll ). NGC 424 has a mid- 
IR spectral index of a —1.1, making it the blues t type 
2 compared to the other objects in ,H6nig e t aD ()2010l ) 
with only one type 1 AGN showing a bluer mid-IR spec- 
trum. This points toward a lack of cool dust relative to 
the amount of hot dust in NGC 424. 

3.2. Interferometry data analysis 

The visibilities and correlated fluxes of all reduced 37 
datasets are shown in the online-only Fig. [101 These data 
will be analyzed in the following. 

3.2.1. The size of the mid-IR source m NGC 4^4 

A note on interferometric sizes — The correlated fluxes 
can be compared to the total fluxes to obtain visibilities 
that, in turn, can be converted into source sizes using 
a suitable model. AGN sources have often been mod- 
eled by means of a Gaussian flux distribution resulting 
in sizes based on the Gaussian FWHM or HWHM. How- 
ever, it was shown recently that sizes inferred from single- 
Gaussian models can be quite misleading since the under- 
lying brightness distribution barely resembles a Gaussian 
(Kishimoto et al. 2011b). In fact, these sizes can be con- 
sidered rather a measure of the array resolution than the 
source size. This becomes even more of an issue when 
2-dimensional sizes are determined from data where the 
baseline lengths change with position angle - meaning 
the source shape is depending on the u?;-coverage, a typ- 
ical problem in contemporary IR interferometry owing to 
the limiting number of telescopes in the arrays. There 
are two ways to approximately overcome this problem. 
First, one chooses a fixed visibility, infer the baseline 
lengths/spatial wavelengths for all position angles at this 
visibility, and calculates a corresponding size. This has 
been done in Kishimoto et al. (2011b) in terms of the 
half-visibility radius Rvo.5 for a fixed visibility of 0.5 or 
related half-light radius Ri/2- Second, it is also possible 
to take a fixed projected baseline length (PEL) and infer 
the visibilities for all position angles at this PEL. The 
first method is more related to the intrinsic structure of 
the object while the second one reflects the shape that 
a telescope with an aperture the size of the PEL would 
see. Eoth methods, however, overcome the problem of 
inhomogeneous wii-coverage to some extent. 

NGC 424 emission source — In the left panel of Fig. [SJ 
we show the position-angle-dependent Gaussian HWHM 
sizes of the mid-IR emission source in NGC 424 at a 
fixed PEL of 58 m as obtained from our data. Where 
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Figure 5. Sizes of tlie mid-IR emission of NGC 424 for different wavelengths from 8.4 fira (red circles) to 12.5 ^m (blue circles), determined 
from our MIDI interferometry data using different methods. Left: Gaussian HWHM sizes (in mas) at a fixed baseline of 58it3 m for a 
range of position angles (PAs). Center: Gaussian HWHM sizes (in mas) at a fixed baseline of lOOitSm for a range of PAs. Where no 
individual data were available, data pairs of similar PA have been interpolated or extrapolated. The colored-dashed lines are elliptical fits 
to the data at the corresponding wavelength bin. Right: Half-light radius for a range of PAs. The black-dashed line is an elliptical fit to 
the wavelength-averaged data. 



Table 3 

Geometric fits to the visibilities of NGC 424. 



fit type 


wavelength 


semi 


major axis 


a 


semi 


minor axis 


b 


position angle 




(fim) 


(mas) 


(pc) 


(nn) 


(mas) 


(pc) 


(nn) 


(deg) 


Gaussian 


9.1 


8.0 ± 1.0 


1.7±0.2 


26 ±3 


6.5 ±0.7 


1.4 ±0.2 


21 ± 2 


-40 ± 23 


(PEL 58 m) 


11.0 


8.6 ± 1.2 


1.9 ±0.3 


28 ±3 


7.9 ±1.0 


1.7 ±0.2 


25 ± 3 






12.5 


10.2 ± 1.7 


2.2 ±0.4 


33 ±5 


8.9 ±1.8 


1.9 ± 0.4 


29 ± 5 






8.4-12.5 


8.7 ±0.9 


1.9 ±0.2 


30 ±3 


7.6 ±1.0 


1.7 ±0.2 


26 ± 3 


-40 ±8 


Gaussian 


9.1 


5.7 ± 1.5 


1.2 ±0.3 


18 ±4 


4.0 ±0.3 


0.9 ± 0.1 


13 ± 1 


-12 ± 26 


(PEL 100 m) 


11.0 


6.5 ± 1.6 


1.4 ±0.4 


21 ±5 


4.9 ±0.6 


1.1 ± 0.1 


16 ± 1 






12.5 


6.8 ± 1.1 


1.5 ±0.2 


22 ±3 


5.7± 1.3 


1.2 ± 0.3 


18 ±4 






8.4-12.5 


6.2 ±0.8 


1.4 ±0.2 


21 ±2 


4.7 ±0.7 


1.0 ±0.2 


16 ± 2 


-5± 11 


Half-light radius 


9.1-12.5 


9.1 ±0.9 


2.0 ±0.2 


31 ±3 


7.0 ±1.1 


1.5 ± 0.3 


24 ±4 


-21 ±3 



no data were available in the range of 55-61 m, we used 
pairs of data at similar position angles (PAs) and in- 
terpolated or extrapola ted to 58m (including errors). 
(|Kishimoto et al.ll2011bO showed that the brightness dis- 
tributions in AGN follow simple power laws, so that we 
linearly inter-/extrapolated in log-space. Because of the 
resulting well-sampled PA coverage, some of the source 
characteristics can be easily inferred. First, the object 
is smallest at the shortest wavelength bin and becomes 
larger with wavelength (see also Fig. [SI). This is qualita- 
tively consistent with the idea that hotter dust (primar- 
ily emitting shorter wavelengths) is located closer to the 
AGN. Second, the mid-IR source is non-circular. While 
this is not too obvious in the 58 m data, we also show 
Gaussian HWHM sizes at 100 ni in the middle panel of 
Fig. [S] The sizes in East/West direction are clearly 
smaller than those close to North/South for all wave- 
lengths. It should be noted, however, that there are 
larger gaps in the PA coverage at 100 m because of the 
limited uw-coverage by the UT telescopes (see Fig.[T]). In 
the right panel of Fig. \5\ we show the PA-dependence of 
the half-light radius in NGC 424. As previously noted, 
this measure is independent of the actual baseline length. 
However, data interpolation is reversed with respect to 
the fixed-baseline method and results in larger uncer- 



tainties of the spatial wavelength because the change in 
visibilities with baselines is small. In the right panel of 
Fig. \5\ we show all data pairs that allow us to establish 
a clear half-light radius. The problem is particularly rel- 
evant to the shorter wavelength data as can be seen by 
the sometimes erratic size changes at 9.1 fim (orange cir- 
cles). Taking all three panels together makes it clearer 
that the objects displays some degree of elongation. 

To be more quantitative about the object's size and 
shape, we fitted ellipses to the different measurements. 
In the left and middle panel of Fig. [5] we show fits to 
each wavelength bin of the Gaussian sizes at 58 m and 
100 m. Because of the higher uncertainties in the deter- 
mination of the half-light radius, only the median visi- 
bility of all wavelengths was used for the fit in the right 
panel of Fig. [5] In Table |3] some reference fitting results 
are listed. The size is typically between 0.9 pc and 2.2 pc, 
depending on baseline length, wavelength, and orienta- 
tion. The half-light size of (2.0 ± 0.3) pc x (1.5 ± 0.3) pc 
would make a variability origin of the difference between 
the MIDI and Spitzer-|-ISAAC fiuxes physically possi- 
ble (see Sect. 13.1]) . but we note again that the difference 
spectrum is cool. Therefore, a final conclusion about the 
origin of the flux difference cannot be made. For the 
sizes please note that due to the itw-coverage, the largest 
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Figure 6. Wavelength-dependence of the Gaussian HWHM size 
along the major axis (filled circles) and minor axis (open circles) at 
a fixed baseline length of 58 m. Please note that the error bars are 
larger than the bin-to-bin variance because they are dominated by 
systematic errors in the visibility. The dashed and dotted lines are 
linear fits to the data at > 9 fim for the semi-major and semi-minor 
axis, respectively. 
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Figure 7. Position-angle dependence of the mid-IR brightness 
distribution, parameterized by a surface density power law with 
index a and a fixed temperature law (see text for details). The filled 
circles mark PAs where the brightness distribution was modeled 
by two or three sets of interferometric observations at different 
baseline lengths, while the open circles denote PAs with just one 
baseline. 



(relative) uncertainties in the sizes are toward the semi- 
minor axis for the 58 m fit and toward semi-major axis 
for the 100 m fit. 

We present the sizes in observed (milliarcseconds) and 
physical (pc) units, as well as in intrinsic units, i.e. in 
units of the sublimation radius. The latter one is very 
useful when we want to compare pro perties in differ- 
ent ob j ects wi th different luminosities (|Kishimoto et al.l 
I2009U I2011bl ). For that we have to establish an in- 
ner radius, r-m of the dust distribution (also referred 
to as sublimation radius). Usually this is inferred ei- 
ther from reverberation rnappiiig or near-IR interferom- 
etry (e .g. iSuganuma et al.ll2006tTKishimoto et al.ll2009al 
1201 lal) . Since we do not have either information in this 
type 2 AGN, we use an e stimate based on t he 12 fj,m 
luminosity as discussed in IHonig et al.l (j2010t) . Taking 
the mean between Spitzer and MIDI total flux, we ob- 
tain Tin ^ 0.065 pc for a luminosity of i^i,y(12/im) — 
5.4 X 10^3 ej.g/s^ 

The ellipse fits illustrate some key characteristics of 
the mid-IR emission source. They confirm the obvious 
trend that the source size increases with wavelength, as 
illustrated in Fig. |6l Although the error bars are large 
(dominated by systematic errors in the visibility), a sys- 
tematic trend of larger mid-IR emission sources for longer 
wavelengths is seen. 

3.2.2. Position angle dependencies of the mid-IR source 

The elliptical fits also provide information of the di- 
rection of elongation. In Table [3] we list the position 
angle of the half-light radius for those fits that allow 
for such a constraint. While the formal fitting errors 
are considerable, all position angles systematically indi- 



cate an elongation toward the NW/SE direction. The 
fits that average over all wavelengths give approximately 
-5° < PA < -40°. There may be a trend that the 
position angle changes slightly from the lower-resolution 
data at 58 m to the 100 m data. While such a trend 
wou ld be reminiscent of i nterferometry data of the Circi- 
nus ([Tristram et al.ll2012f ). the uncertainty in the geomet- 
ric determination of the PA do not allow us to formally 
establish such a skew (see also the PA for the half-light 
radius), in particular because the 100m data do not have 
any coverage in the direction of the major axis. There- 
fore, any position angle determined for the major axis in 
the 100 m data must be considered uncertain. 

We note that this elongation is not affected by the in- 
homogeneous MZ)-coverage, since this has been taken care 
of as discussed above. However, when deriving sizes for 
each wavelength, only the spatial information of the in- 
terferometric data is used. We will now include the spec- 
tral information in our analysis to test the result from the 
geometric model fitting. For that we use a simple model 
of dust emission to characterize the mid-IR emission at 
as many position angles as possible that are covered by 
our observations. For each PA, we simulate a brightness 
distribution of the form 

5.(r) = S.(r(r))./o(^^y (1) 

where r is the distance from the AGN, rin the inner ra- 
dius of the dust distribution, /o is the normalized sur- 
face density at the inner radius, a is the surface density 
power law index, and the black-body radiation with 
temperature T(r) at distance r. In the following analysis 
we will fix /o = 1 for simplicity but use it in Sect. |4]as 
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a free parameter. For the temperature we assume that 
the dust is in local thermal equilibrium and parameter- 
ize the temperature gradient by T{r) = Tguh ■ (''/Hn)^- 
We use a sublimation temperature Tguh = 1 500 K and 
ad opt = —0-36 fr om the near-IR dust opacity model 
bv iBarvainisI (|1987[ ). The only remaining parameter is 
the powe r-law index a given the r^in value estimated in 
Sect. 13. II This optically-thin model describes the directly 
AGN-heated dust as we see it in the mid-IR, while any 
possible cooler dust is not covered. However, the far-IR 
SED does not suggest the presence of such a cool com- 
ponent (see Fig. [Hand Sect. 13. 

In Fig. [7] we show the PA-dependence of the surface 
density distribution power law index a. Where applica- 
ble we fitted the model to two or three data sets with dif- 
ferent baseline length (filled circles) . If no second set was 
available within 10° single baseline data was used (open 
circles). The data clearly show a change of a with PA. 
From about PA —90° to —45°, the power law index be- 
comes less negative. Toward larger positive PAs, a turns 
more negative with a peak in the range of 50 — 80°. This 
change can be either interpreted as a PA-dependence of 
the brightness distribution (i.e. traced back to a change 
in a, /3, or a combination of both), or a PA-dependence 
of Tin, which is used for normalization. In either case it 
illustrates that, indeed, the mid-IR source changes struc- 
turally with position angle, and the shallowest a values 
(=least negative) in PA —30° corresponds to about 
the same direction as the major axis in the geometric 
fits. 

4. SIMULTANEOUSLY MODELING THE SED AND 
INTERFEROMETRY OF NGC 424 

4.1. Model description 

The analysis of the NGC 424 interferometric data 
presented so far was based on the observations with 
as little model-dependence as possible in order to ex- 
tract basic properties of the mid-IR emission. How- 
ever, only specific aspects of the data have been investi- 
gated. Now, we want to model the full data set of mid- 
IR interferometric observations together with the SED. 
This can be done, in principle, with advanced models 
of clumpy tori (for such exampl es s ee iHonig et al.l [20061 
120101: iTristram fc"Sc hartmann 20Tl|) . On the other hand, 
the visibilities in NGC 424 do not go below 0.4 and mod- 
els suggest that effects that clearly distinguish homoge- 
neous models from clumpy on es will only appear when 
an object is better reso lved (jSchartmann et all 120081 : 
IHonig fc Kishimotol IMol) . This means the present ob- 
servations are not sensitive to substructure but only to 
the overall distribution of the brightness. 

We, therefore, decided to use a rather simple physical 
model to reproduce the interferometric and photometric 
data simulta neously to avoid over-int erpretation. The 
model follows iKishimoto et al.l (|2011b[ ) and was success- 
fully applied to type 1 AGN. The brightness distribution 
is defined by the dust surface density power law index 
a for the emitting dust and the temperature gradient /3 
of the dust (see Eq. [T|), plus an additional hot ring com- 
ponent (see below). Here, [3 is left as a free parameter 
and T(r) starts at a maximum temperature Tmax that 
we also fit for. The surface density fo at the inner radius 
will be used as a free parameter now (see Sect. I3.2.2p . 



We note that we do not need an ad-hoc assumption of 
Tin but can leave it as a free parameter as well. It sets 
the normalization of both the fluxes and the visibilities: 
For a given set of parameters, a larger r^^ will result 
in higher fiuxes (larger emitting surface) and lower vis- 
ibilities (larger sizes). As already discussed, the mid-IR 
emission is elongated. Therefore, we describe the source 
shape, and as a consequence the brightness profile, as an 
ellipse with the ratio s — a/b between the semi-major 
and semi-minor axis and orientation 4> of the major axis 
as parameter^ 

The hottest emission in this model originally designed 
for type 1 AGN is coming from a ring (in addition to 
the power law) at about the inner radius with tem- 
perature Tin = 1 400 K and with a contribution fac- 
tor /hot- This hot ring has been held responsible for 
the observed near-IR b ump in the SEDs of many type 
1 AGN (e.g. lEdelson fc MalkanI [19861: IKishimoto et al.l 
I2011bl: IMor fc Netzedl201lD . Since NGC 424 is a type 2 
AGN, we will show model fits with and without the hot 
ring. The origin of this hot emission component is not 
fully understood. It may either be the result of a differ- 
ent dust size/composition in the hot region of the torus 
(silicates sublimate at cooler temperatures than graphite 
grains; small grains heat up to the sublimation temper- 
ature at larger distances), or coming from a graphite- 
dust component in the innermost part of the polar out- 
flow region. We will discuss an addition possibility in 
Sect. 15.41 In any case the near-IR bump is typically not 
seen in type 2 AGN, indicating significant extinction to- 
wards this component. From the presence of the silicate 
features in emission, we conclude that the dust column 
toward a hot dust component may be optically thin at 

> 8 /im. From about 4—8 /im typical ISM dust extinction 
curves do not change significantly (but note that the dust 
may have a non-ISM composition; see Sect. 15.4]) . How- 
ever, towards the near-IR J, H, and K bands - the peak 
regions of the hot dust component - the extinction in- 
creases dramatically by a factor of 5 — 10. To avoid these 
subtle problems, we ignore the shortest wavelengths in 
the fitting process. 

4.2. Data selection 

As discussed in the previous section, there may be sev- 
eral problems involved when trying to fit the near-IR 
photometric data, most notable the possible extinction. 
We, therefore, restrict our modeling efforts to the data 

> 3 fiiR. In order to capture the downturn of the SED 
longward of 20 /zm, but not overweighting the photomet- 
ric data over the interferometric, we extract fluxes at 
22 fim and 32 fim from the Spitzer data and treat them 
as photometric data points. In addition, to account for 
the offset of the MIDI data with respect to the rest of the 
SED (see Sect. 13. 1|) the photometry has been corrected 
by a factor of 1.3, so that the Spitzer IRS 8 — 13 /im 
matches the MIDI spectrum. 

To speed up the fitting process and gain SNR, we 
resampled and binned the interferometric data. The 

Since NGC 424 is a type 2 AGN, probably at considerable 
inclination, the meaning of a ring as description for a projected 
face-on torus as in a type 1 AGN becomes obsolete. The concept 
of an inner radius with inclusion of ellipticity, however, captures 
the fundamentals of dust emission and sublimation and axial sym- 
metry. 
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Table 4 

Best-fit parameters for the modeling of the IR photometry and interferometry. 

model a P 4> s = a/b T„,ax /o 7hZ ~ xl 

(dcg) (K) (pc) 

photometry -|- interferometry 

hot ring + T/d-gradicnt -0.97 ± 0.05 0.31 ± 0.03 -27 ± 2 2.2 ± 0.1 570 ± 63 0.81 ± 0.20 0.05 ± 0.02 0.101 ± 0.013 0.86 

T/d-gradicnt (w/o L, Af) -0.85 ± 0.14 0.32 ± 0.03 -27 ± 1 2.2 ± 0.1 662 ± 81 0.99 ± 0.25 ... 0.066 ± 0.013 0.83 

T/rf-gradicnt (MIDI only) -1.00 ± 0.18 0.31 ± 0.04 -27 ± 2 2.2 ±0.1 713 ± 101 0.94 ± 0.30 0.063 ± 0.016 0.80 

photometry only 

hot ring + T/d-gradiont > -0.5 < 0.36 ... ... > 900 < 0.3 >0.1 < 0.05 

T/d-gradicnt (w/o A/) -1.02 ± 0.25 0.32 ± 0.08 ... ... 545 ± 154 unconstr. ... unconstr. 0.32 

Note. — Model parameters: a . . .power law index of the surface density distribution; /3 . . .temperature-gradient power law index; (p . . .position 
angle of the major axis (in degrees from N toward E); s . . .ratio of major axis a and minor axis b; Tmax . . .maximum temperature of the dust 
distribution; /o . . .normalized surface density at the inner radius rin; /hot . . .fractional contribution of hot-dust component; rin . . .inner radius of 
the dust distribution 



8 — 13 /iiTL range will be represented by 9 wavelength bins 
per data set in the fitting process (see Fig. [10] in the 
appendix). While the 2009 and 2010 data are in gen- 
eral agreement for all baseline configurations that match 
between both epochs, two of the 2007 data sets show 
significant disagreement. In particular, the two U3-U4 
sets with only moderately different PBL and PA (see 
Table. [T]) are inconsistent with each other. There are 
several technical difficulties involved when dealing with 
the 2007 data (e.g. the IRIS beam-tracking system of 
the VLTI did not work properly at that time on tar- 
gets with X-magnitude ^ 10, resulting in poor/unstable 
beam overlap) in addition to the low-quality total flux 
reference of 2007 (see Sect. 12. ip . so that we will not in- 
clude these data in the fitting process. However, the fits 
will be overplotted onto the 2007 data in Fig. [TU] as a 
comparison. 

Last, the visibilities and correlated fiuxes show a dis- 
tinct change of slope at the shortest wavelength bin at 
about 8.4 /im with respect to the rest of the spectrum. 
There are two possibilities that might explain this feature 
in the data. First, the change of slope could be real. As 
a consequence the wavelength-dependent sizes would not 
become smaller with decreasing wavelength anymore. In- 
deed the 8.4 /im sizes may be even systematically larger 
than for the neighboring 9.1 /im bin (see Fig. IH). While 
this is difficult to explain with in simple radiative equi- 
librium, it could be an effect of extinction/obscuration at 
shorter wavelengths. However, as discussed in the pre- 
vious section, the overall increase of extinction toward 
shorter wavelengths in this range is rather small for any 
kind of dust (ISM, large grains, silicates, graphites). An- 
other possibility for a real decrease in visibility is con- 
tamination by the 7.7/im PAH feature. Indeed, this 
feature is weakly present in the Spitzer IRS spectrum, 
and the wings may still be seen on the short-wavelength 
side of our MIDI total flux spectrum. If these wings are 
not present in the correlated fluxes, but only in the to- 
tal flux, the visibility would be expected to be lower, 
as seen here. However, in this case the lower visibilities 
are not physically related to the parsec-scale dust distri- 
bution. Second, the lower 8.4 /im visibilities/correlated 
fluxes may be an instrumental/calibration effect. Indeed 
the sensitivity of the MIDI detector is sharply changing 
at these short wavelengths, in addition to generally lower 
flux levels than at longer wavelengths and shorter coher- 
ence times. To avoid problems with these lower visibility 



at 8.4 /im that are most likely not related to the parsec- 
scale dust, we will exclude it from the fitting process (i.e. 
using exclusively 9 — 13 /im data). 

4.3. Model results 

In Table |4] we list the results for a number of different 
modeling approaches (models as described in Sect. 14.1] ). 
The first model uses all the data as described in Sect. 14.21 
and fits the temperature/density gradient model plus the 
hot ring. The second fit excludes the L- and Af-band 
photometry as well as the hot ring that is essentially only 
necessary to account for these wavelengths (see Sect. 14. 1( ). 
The resulting reduced Xr of both models are equal, un- 
derlining this point. Finally we use only the MIDI data in 
the 9 — 13 /im range and model it with the temperature- 
and density- gradient (or T/d-gradient) model without 
hot ring, again with a Xr ^ 0-8 as the previous fits. 

The resulting parameters point toward a surface den- 
sity distribution with power-law index a — 0.9 ... — 
1.0. This result is consistent with the more model- 
independent determination of a illustrated in Fig. [71 tak- 
ing into account that the more physical model accounts 
for the PA-dependence by ellipticity. The value of a is 
rather on the steeper sid e when compared to type 1 AGN 
(jKishimoto et al.ll2011b[ ). which is particularly interest- 
ing since torus models would predict generally redder 
colors for type 2s than type Is because of obscuration 
effects. 

The orientation of the mid-IR emission is quite sim- 
ilar to what we got from the geometric fits (consistent 
within error bars). A PA of -27° ± 2° is within the 
range of PAs as listed in Table [31 The a/b ratio of 2.2 in 
the model illustrates the strong ellipticity of the mid-IR 
source. This ratio may seem larger than visible in Fig.[5l 
but these perceived differences come, at least in part, 
from resolution effects (see also the stronger ellipticity 
for the 100 m data than for the 58 m datafl- 

For illustration we also fitted the model to the 
photometry-only data. The model using the hot ring 
+ temperature gradient is basically unconstrained with 

^ Similar to single-telescope data, the observed/measured object 
shape depends on the intrinsic shape convolved with a "resolution 
kernel" or PSF. The size of this kernel depends on the baseline 
lengths and has larger widths for smaller baseline lengths. This 
leads to "smearing" of the actual shape of the object. In an un- 
resolved case, the observed shape would be the PSF of a point 
source. 
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Figure 8. 11 /jm visibilities plotted against the projected baseline 
length for all 2009/2010 observations except the UT1-UT3 data 
(see note in the text). The position angles for each data point 
is color coded (from dark-blue along the modeled major axis at 
—27" to red along the minor axis in PA 63"). Overplotted are 
PA-dependent visibilities for the T/d-gradient model (w/o L/M 



photometry) as described in Table. 



the range of parameters of 2 x Xr-min given in Table ID 
The same happens when using the T/d-gradient model 
on the MIDI-only photometry (i.e. the 8 — 13 /im to- 
tal flux spectrum). Both limits are, at best, marginally 
consistent with the fits involving interferometry. Bet- 
ter constraints are achieved for the T/d-gradient model 
fitted to the mid-IR photometry without the L- and M- 
band fluxes. Although the range of parameters allowed is 
still considerable, it is in overall agreement with the fits 
including interferometry. This illustrates that the spa- 
tial information is important to properly constrain such 
models. Using complex models on SED-only data may, 
therefore, lead to results contradictory to the actual spa- 
tial distribution of the dust, unless the wavelength range 
and/or parameters are carefully chosen (for a more de- 
tailed discussion with resp ect to clumpy torus models see 
iHonig fc KishimotollMO) . 

Finally, we show the 1 1 fim visibilities versus baseline 
in Fig.|8]for all baseline configurations, except UT1-UT3. 
The data sets for this configuration have been taken un- 
der the worst observing conditions, which resulted in dif- 
ficult calibration, large error bars, and the lowest SNR 
of all configurations (see Table. [T]). Therefore, we ex- 
clude these sets from the plot for clarity. The PA for 
each data point is color coded to be able to distinguish 
the baseline configurations closer to the major axis (to- 
wards dark-blue) and to the minor axis (towards red) as 
set by the model orientation. As discussed before, the 
inhomogeneous distribution of telescope position angles 
and projected baselines make a direct assessment of any 
kind of elongation difficult. However when concentrating 
on the shorter baseline lengths (50-60 m), a trend is seen 



that the blue data points show systematically lower vis- 
ibilities than the green points. At the longest baselines 
(>100m), the green data points, about in between major 
and minor axis, are also systematically lower than the red 
data that are pointing approximately in the direction of 
the minor axis in the model. As a comparison and further 
way to prove the elongation and direction of the mid-IR 
source, we overplotted the visibility curves for different 
position angles based on the T/d-gradient model without 
L/M-band photometry. This illustrates that, in spite of 
the scatter and the inhomogeneity of the data, the mid- 
IR source shape and radial brightness distribution are 
well reproduced by the elliptical T/d-gradient model. 

5. DISCUSSION 
5.1. Mid-IR emission elongated in polar region 

The orientation of the major axis can be put in the 
context of the AG N system axis in NGC 424. Although 
iNagar et all ()1999D report that their 20 cm/ 1.5 GHz radio 
observations may indicate a marginal extension in about 
the same direction as their beam size, this could not be 
confir med with their higher resolution 3.6 cm/8.4 GHz 
data. iMundell et al.l (|2000( ) consider the source to be 
unresolved at both frequencies in t heir data with bette r 
SNR, in fine with the analysis bv iThean et al.l ()2000D . 
Consequently, a possible linear radio structure that could 
i ndicate the system a x is is n ot detected. 

iDurret &: BergeronI (jl987[ ) report [O III] imaging ob- 
servations of the nucleus. They note an elongation of 
the [O III] in PA^ 70°. The extension is seen best in 
the faintest two contours of 1.25 and 2.5% relative to 
their peak contour. The authors note that the determi- 
nation of the continuum (dominated by the host galaxy) 
is uncertain, which may result in problems subtracting 
the continuum emission from the line image. They give 
an elongation of the [O III] emission in PA 70°, with the 
galaxy major axes in PA 59° (i.e. APA = 11°). Unfor- 
tunat ely, no HST im age at higher quality is available. 

Mo ran et al.l (|200C1( ) report results from optical spectro- 
polarimetry of the nuclear region of NGC 424. While the 
unpolarized spectrum shows signs of a strong host com- 
ponent of 45% (quite typical for Seyfert 2 AGN), the po- 
larized flux displays broad emission lines, which indicates 
a hidden type 1 nucleus. From their polarization data, 
the authors give a polarization angle of 43°. In type 2 
AGN, the polarization angle has been found to be domi- 
nated by scattering (either electrons or dust) in the polar 
region of the AGN and be approximately orthogonal to 
the l inear radio axis where such a feature was dete cted 
(e.g. IAntonuccil[l98l [l993l [200l ISmTth et al.ll2004]) . It 
is, therefore, also orthogonal to the symmetry axis of the 
torus. In order to confirm that the scatterer is located 
in the polar region, we es timated the intri i isic p olariza- 
tion of the scattered fight. iMuravami et al.l ()1998[ ) report 
deep spectral observations of the Ha line and detect shal- 
low broad wings to the otherwise narrow emission line. 
We extracted the flux from their spectru m and co mpared 
it to the polarized (Stokes) flux from . Moran et al.l ([SOOO). 
The broad component of the Ha line shows about 15% 
polarization, which is typical for scattering over a siz- 
able solid angle as in the polar region of the AGN (e.g. 
lAnto nucci fc Milleiill985l ). 

Based on the polarization measurements, the system 
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axis for NGC 424 would be established in approximately 
PA 43° - 90° = -47°. This is not consistent with the 
[O III] PA. Both methods are not a ccurate, though 
and offsets have to be expected (e.g . lAntonuccl 120011 : 
iSchmitt et alll2003l : ISmith et al.ll2004[ i. In the particular 
case of NGC 424 - based on our analysis of the polar- 
ization data and given the mentioned uncertainties with 
the [O III] data - we consider the polarization data as the 
more reliable method to estimate the system axis. As a 
consequence, we expect that the system axis is close to 
PA —47°, which is very close to the direction of the ma- 
jor axis of the mid-IR emission (PA —27°), meaning that 
the parsec-scale mid-IR emission shows elongation in the 
polar direction. 

5.2. Polar elongation in other AGN 

The elongation of the mid-IR emission in the polar 
direction, instead of the "disk/torus plane" , has been 
observed recently in some objects fe.g. iBock et~alll2000l: 
iPackham et"alll200l . 

Circinus galaxy: — IPackham et al.l (|2005l ) shows mid-IR 
images of the AGN in the Circinus galaxy at 8.8 //m and 
18.3 /im with some extended emission ~40pc to the East 
and West of the nucleus that coincides s patially with the 
rim of the known [O III] outflow cone. iReunanen et al.l 
([2OT0h compare their VLT/VISIR 11.8 and 18.7yLtm im- 
ages to near-IR [S VII] coronal line images and demon- 
strate that the mid-IR extended emission traces the NLR 
region. fPackham ct al. (2005) further note that this ex- 
tended emission is responsible for about 30% of the total 
emission in the central 3'.'5 while the other ~70% orig- 
inate from the unresolved nucleus (<20pc in their im- 
ages). Both the color temperature and intensity of the 
extended emission are in good agreement with optically- 
thin dust in local thermal equilibrium and heated by the 
AGN. 

NGC 1068:— The outflow region of NGC 1068 dis- 
plays significant mid-IR emission tha t extends out to 
more than 50 pc fro m the nucleus (e.g. iBock et 311120001 : 
iTomono et al.ll2001f) and associated w ith emission from 
dust clouds fe.g. iTomono "eFaLl [20061 ) ■ The extended 
emission contributes about 70% to th e 11.6 ixm flux in 1" 2 
or 85 pc aperture around the nucleus (iMason "eraLl[2006h . 
Yet, on parsec scales the emission is dominated by the cir- 
cumnuclear dust and exte nds about perpendicular to the 
large -scale outflow cavity (jJaffe et al.ll2004l : iRaban et al.l 
l2009f ). It is noteworthy, though, that the elongation in 
outflow direction has also been observed by Speckle imag- 
ing in the near-IR [H and K ba nds) with a iC-band size 
of l.Spc X 2.8pc (e rror ±0.3pc: IWittkowski et 311119981 : 
IWeigelt et aLl[200l . While one m3y 3rgue 3bout pecu- 
li3rities of NGC 1068, NGC 424 3nd Circinus indic3te 
thst a polar-extended mid-IR emission is not untypical 
in AGN. 

Other objects: — iHonig et al.l (|2010r ) presented mid- 
IR 8m-telescope images of IC 5063 and MCG-3-34- 
64 that also showed mid-IR emission in polar direction. 
In these cases, however, the extension is observed on 
scsles of >100pc, while in NGC 424 we are dealing 
with a factor o f 100 s maller dist3nce. In NGC 1386, 
lReun3nen et al.l (|2010[ ) report extended mid-IR emis- 
sion in approximately North/South direction. This is 



at 3bout the same direction as the [O III] outflo'v^O. 
( Ferruit et all I2000D and radio -emission position angle 
( Nag3r et 3l.lll999i rin this object. Finally, we note that 
the type 1 AGN NGC 3783 is another candidate with 
even s tronger elongation in polar direction on parsec 
scales (jHonig et 3l . ,2010t iKishimoto et 3l.ll2011bl Honig 
et 3I. 2012, in preparation). 

5.3. Where is the torus? 

The extended mid-IR emission in polar region in 
NGC 424 gives rise to the question if it is consistent 
with the paradigm that emission from the dusty torus 
domin3tes the mid-IR emission. This paradigm is well 
summarized in Mor et al. ( 2009 , their Fig. 4) where the 
3uthors show results from decomposition of Spitzer SEDs 
of type 1 AGN into the torus, a dust NLR, and an ad- 
ditional hot dust components. According to the median 
of all decompositions, the torus contributes >80% of the 
flux between 7 and 15 /im and dominates (>50%) in the 
range from 4 to 25 /im (the long- wavelength side being 
limited by wsvelength coverage of the data). However, 
our data suggests that in NGC 424 dust in the polar 
region dominates the mid-IR emission. Indeed, we do 
see elongation in polar direction at all observed baseline 
lengths, i.e. a consistency in elongation toward polar di- 
rection at these spatial scales. At the longest baseline 
lengths, the visibilities have dropped to ~0.4, meaning 
that at least 60% of the mid-IR single-aperture/total flux 
show this elongation direction. Thus, it is not too iai 
resching to interpret this as evidence that > 60% of the 
mid-IR total flux emission originates from the polar re- 
gion. 

Interestingly, iTristram et al.l (|2012[ ) report on new 
mid-IR interferometric observations of the Circinus 
galaxy on short baseline lengths, com plementin g pre- 
vious results on longer baselines (Trist ram et al.l [2007'). 
These new dat3 indicate that the mid-IR emission is 
extended in outflow direction not only in the single- 
telescope aperture as mentioned above, but 3lso on the 
new bsselines. Only 3t baseline lengths longer than 
~50 m does the orientation of the mid-IR emission (de- 
scribed as 3 0.4 pc "disk") chsnge direction toward per- 
pendicular to the outflow region and, therefore, in the 
direction expected for the classical torus. However, we 
note that at these baselines the visibilities decreased to 
'^0.15, meaning that 85% of the single-aperture mid-IR 
flux originates from the pol3r region on sc3les 3S smsll 
3S few parsecs. 

Whether our dat3 for NGC 424 can be considered con- 
sistent with the paradigm of mid-IR emission primarily 
originating from the torus depends on the question if 
torus models can reproduce the elong3tion in polar di- 
rection in an obscured AGN. It might be possible thst 
self-shielding effects in the mid-pl3ne of the torus C3use 
the torus to appear elongated in polar direction while the 
dust is still strongly concentr3ted to the mid-plane. An 
addition3l f3ctor could be 3nisotropic illuminstion from 
the BBB/accretion disk onto the torus, i.e. radistion 

* We loosely refer to the narrow emission line regions 
(NLR) as outflows based on the evidence from both statisti- 
cal analysis of emission lines in object samples and individua l 
sources (e.g. Storchi-Bergmann ct al. 1992; IReunanen et al1l2003l : 
.Miill cr-Sanchcz ct al. 2011.; Antonucci 200i, Fig. 7) 
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escaping preferentially in polar direction. Detailed ra- 
diative transfer models that also include visibility studies 
have been presented by 'Schartmann et al.' (2008, smooth 
and clumpy) and Honig & Kishimoto (2010, clumpy), 
the former also invoking anisotropy of the accretion disk. 
Both studies agree that for sets of models that do well 
represent the IR SEDs as observed for typical AGN, 
the visibilities should be higher in polar direction (i.e. 
smaller size) than along the mid-plane of the torus in 
type-2-like edge-on viewing angles, and more circular to- 
wards face-on orientations. Strong polar extension is not 
to be expected. This also means that although there 
is self-shielding in the equatorial plane of the torus, the 
dust mass in typical radiative transfer models is not suffi- 
cient to cause strong-enough obscuration that the mid-IR 
emission is seen elongated in polar direction. 

The radiative transfer models use ad-hoc parameter- 
izations of the dust distribution, while hydrodynamical 
simulations attempt at getting the mass distribution self- 
consistently fro m global parame ters of the AGN and its 
environment. iWa da et al.l ()2009f) present high-resolution 
numerical simulations of ISM accretion in the inner 32 pc 
of an AGN. Although the authors did not perform radia- 
tive transfer simulations of the hydrodynamic models, 
the mass density distribution in the inner parsecs is quite 
flat (disk-like). At larger radii (>10pc), the disk flares 
resulting in a scale hight of the order unity, which is 
also typically used in the previously mentioned radiative 
transfer models. Therefore, it is not expected that these 
models result in strong p olar-extended mid-IR emission. 
iSchartmann et all (|2009[) present hydrodynamic simula- 
tions of the central 100 pc around an AGN, assuming 
that the dust is produced in-situ by a stellar cluster. 
The authors show (log-scale) mid-IR radiative transfer 
images and SEDs of their standard model. Indeed, these 
models do show elongation in polar direction and might 
be considered a better match to our observations. As 
the authors note, the polar elongation is the result of 
heavy self- absorbing in a dense torus-plane region, so 
that only the funnels of their mass distribution are il- 
luminated. As a consequence, the overall IR emission 
is dominated by cool dust resulting in very red SEDs, 
incompatible with NGC 424 (note that the SED of the 
type 2 AGN NGC 424 is significantly bluer in the mid- 
IR than those of type 1 AGNs NGC 4151 and Mark 841 
that have been modeled using_^these simulations; see Fig. 
10 in Schartmann et al. 2003). In particular, the ob- 
served silicate emission feature in NGC 424 would require 
a orientation close to face-on in these models, which, in 
turn, would make the mid-IR emission appear rather cir- 
cular (i.e. dominated by the inner disky structure as 
noted by the authors) and probably inconsistent with 
the a/h = 2.2 axis ratio in NGC 424. Additionally, the 
observed high Hydrogen column density would be diffi- 
cult to reproduce. 

As mentioned above, self-absorption of mid-IR radia- 
tion in the plane of the torus might be a possibility for 
the polar extension of the mid-IR emission. If the polar- 
extended mid-IR emission was caused by a high optical 
depth in the mid-IR (i.e. the mid-IR emission of the 
dust in the torus gets reabsorbed), then we would have 
to expect that the "missing" mid-IR emission from the 
plane is re-emitted at longer wavelengths. Consequently 
we would see either a small bump in the nuclear SED 



at longer wavelengths (where the optical depth becomes 
smaller than unity), or, at least, that the nuclear SED 
appears very red (as illustrated by th e hydrodynamic 
simulations of ISchartmann et al.l [2009( ). In NGC 424, 
however, neither is seen. On the contrary, longward of 
the MIDI-covered 8 — 13 /xm range and the 18 /im silicate 
feature, the SED turns down (see Fig. meaning that 
there is probably no significant amount of radiating dust 
at longer wavelengths outside of what we see with MIDI. 

5.4. A possible change of paradigm 

Concluding our analysis from the previous section, it 
is very unlikely that in NGC 424 the elongation in polar 
direction is (directly) associated with the dusty torus. It 
seems much more likely that >;60% of the mid-IR emis- 
sion is generic to th e polar region. When also consid- 
ering the results by iTristram et al.l (|2012j ) on Circinus, 
we may argue that the torus becomes energetically in- 
significant in the mid-IR in these two sources (and prob- 
ably NGC 3783). Considering that a few other sources 
also show extension in polar direction in the mid-IR on 
larger scales (and, in contrast, there are no clear cases 
with dominant disk- like mid-IR emission) , it is not reach- 
ing too far to discuss possible implications for AGN in 
general. 

First, we will have to address the nature of the emis- 
sion in the polar region in NGC 424. As indicated by 
our model results shown above, the emission is consis- 
tent with dust in thermal equilibrium with the AGN. 
It is, therefore, reasonable to assume that the mid-IR 
emission is indeed originating from dust. Since we know 
that >60% of the mid-IR single-aperture emission orig- 
inates from the polar region, we can use the IR SED in 
this wavelength region as a further constraint. The IRS 
and MIDI spectra show weak silicate emission features 
at 10 and 18 /Ltm, which by its own points toward dust 
as the origin. Moreover, it indicates that the emitting 
dusty medium is optically thin to its own thermal ra- 
diation (see similarities and small anisot ropy of m id-IR 
emission of type 1 and type 2 AGN Ho nig et al.i r2011). 
Therefore, if we would look through the polar region onto 
the AGN (type 1 case), the overall column would still be 
optically thin. This would be consistent with the unifica- 
tion scheme that requires mostly unobscured views of the 
central engine along face-on line-of-sights (at least in a 
statistical sense). However, if the dust is optically thin in 
the mid-IR, the weakness of the silicate features would 
argue for a deficiency of silicate grains with respect to 
graphite grains in ISM dust mixtures, as expec ted when 
dust is exposed to the BBB radiation field (e.g. iPhinnevI 
fl98^. 

NGC 424 is a type 2 AGN and it certainly has an 
obscuring medium along the line-of-sight to the accre- 
tion disk and BLR given the observed high X-ray column 
density in spite of displaying a silicate emission feature 
(a mismatch in column density versus optical/IR obscu- 
ration is common in AGN; Maiolino et al. 2001). This 
means that it offers support for the validity of the unifi- 
cation scheme in this source. The unification scheme also 
requires a geometrically-thick torus - for which we did 
not find evidence in the mid-IR here. However, angle- 

^ We note that the photometric apertures involved in the data 
reduction do not cause any artificial cut-off at these wavelengths. 
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Figure 9. Sketch of the dust distribution that contributes to the 
near- and mid-IR emission around AGN hke NGC 424 and Circi- 
nus. The classical torus that provides the angle-dependent obscu- 
ration is geometrically-thick only in the inner region that radiates 
primarily in the near-IR. In the mid-IR the emission is dominated 
by dusty clouds or filaments in polar region that cover a larger 
surface compared to the torus at these wavelengths. The dashed 
and dotted annuli indicate the approximate scales from which the 
near- and mid-IR emission, respectively, originate. 



dependent obscuration does not require that the scale 
height is large at all distances from the AGN, although 
this is what most torus models assume. It is well possi- 
ble that the geometrically-thick part that gives rise to the 
required torus covering factor in the unification scheme 
is limited to the innermost dusty region. Indeed, as 
mentioned before, type 1 AGN show a distinct near-IR 
"3 — 5 /im bump" in the SED from hot dust, and this 
component has diff erent emiss i on ch a racteristics than 
the mid-IR emission iMor et all (|2009[ ): iKishimoto et al.l 
((2011b). The size of this component is consistent with 
AGN-heated, large graphite grains nea r their sublimation 
temperature (jKishimoto et al.ll2011af ). Moreover, IR in- 
terferometry suggests that this hot-dust component has 
at least a comparable or even l arger surface filling facto r 
as the cool dust (see Fig. 11 of lKishimoto et al]|2011bl ). 
One interpretation of this result is a considerabl e cover- 
ing fa ctor of the hot dust component (see also iMor et al.l 
I2OOI for an SED-based argument on the covering fac- 
tor of the hot component). We may, therefore, speculate 
that the classical, obscuring torus is, in fact, related to 
the hot-dust emitting dust component rather than the 
mid-IR emitting component. 



In Fig. [S] we show a sketch of the dust distribution 
around NGC 424 as this scenario would suggest. The 
torus flares strongly in the innermost part, but becomes 
geometrically thin in the outer regions. In this picture, 
the near-IR emission originates from the obscuring part 
of the torus while the mid-IR emission comes primar- 
ily from the polar region. The IR-emitting sources in 
the polar regions have been drawn as clouds, but they 
could also be filamentary. This inhomogeneous structure 
of the dust is motivated by the mid-IR intcrferometry 
of Circinus (Tristram et al. 2007), but also mid-IR imag- 
ing of NGC 1068 (e.g. Gal hano et al.lfeOOSD . Moreover, 
in Circinus (but not necessarily in NGC 1068), the ex- 
tended mid-IR emission in the single-telescope images 
seem s to be associa ted with the edge of the outflow re- 
gion (jTristram et al . 2007, 2012). We also see an offset of 
20° between mid-IR position angle and the polarization 
angle in NGC 424. This could be a sign of an edge- 
brightening effect as expected in a cone-like region filled 
with optically-thin dust. Therefore, we placed the polar- 
region clouds preferentially around the edges of the open- 
ing cone in the sketch. Note that this is not a full picture 
of the distribution of the dust but only a sketch for the 
IR-emitting material and the polar region may well be 
filled more homogeneously with dust clouds. 

5.5. Further implications 

This picture also has several implication related to 
the physical mechanisms around the AGN. The new ob- 
servations presented here can be interpreted as strong 
support for a dusty wind (as us ed in models for the X- 
ray emission of some AGN: e.g.. iKomossa fc Finklll997l : 
iKomossa fc Badi [1991 iTurner et all l2003| ) . This wind 
could be driven by radiation pressure on dust grains 
in the AGN environment. For an AGN in the Seyfert 
regime radiating at about a tenth of its Eddington lu- 
minosity -/jEdd (defined by electron scattering) , radiation 
pressure is sufficient to unbind dust from the gravita- 
tional potential (the corresponding dust Eddin gton lumi- 
nositv i s ^KHd Hnst, 10-^-'^LKHH:[ Pier fc Kr olik 1992]; 
iHonig fc Beckertl[2007h . The mid-IR spectra also hint 
toward the source of the dust that is driven outward 
in polar direction: as discussed in Sect. 15.41 the dust 
seems to be deficient of silicate grains. Silicates have 
a lower sublirn ation temperature than graphite grains 
(iPhi nnevI Il989f l. making them disappear from a dusty 
medium at generally larger distances from the heating 
source than graphite grains (for a given grain size). This 
is consistent with the near-IR intcrferometry and rever- 
beration mapping sizes in AGN that favor large graphite 
grains for the hottest dust at the sublimation tempera- 
ture (Kishimoto et al. 2007, 2009a), and similar require- 
ments from optical extinction curves in type 1 AGN 
(|Gaskell et al.ll200l . This inner region is also the one 
closest to the AGN where the dust is directly exposed 
(i.e. without significant self-obscuration) to the strongest 
radiation field. It is, therefore, plausible to expect that 
the dust is driven away from this inner part of the torus. 
This scenario seems qualitatively in agreement with re- 
cent radiati on feedback sim ul ations in hyd r odyna mic 
simulations (|Roth et al.ll2012D . Keating et all ()2012D re- 
cently presented IR SEDs of wind models that could 
eventually serve as a starting point for a more quanti- 
tative test of dusty winds. 



14 



Honig et al. 



Radiatively-driven winds or radiation pressure in both 
the optical and IR have also been suggested as a mech- 
anism to supp ort the appar ent geom etrical thickness of 
the torus (e.g. iKonigl fc Ka rtic 1994; IeI^ [2000: Kroli3 
[2003). Of course, this mechanism works best in the inner 
part of the torus close to the sublim ation radius b ecause 
of self-shielding effects further out (jKrolikl [20071) . As a 
consequence the geometrical thickness would be limited 
to just the inner few sublimation radii distances from the 
AGN where, incidentally, the dust is hot enough to cause 
th e near-IR bump s een in type 1 AGN. As pointed out 
in lKishimoto et all ()2011b( ). an analogous near-IR bump 
in the dusty disks of young stellar objects (YSOs) has 
been described as a "puffed- up inner rim" (jNatta et al.l 
[200lHDullemond et al.|[200l . In this scenario we would 
see the near-IR bump pronounced or weak/absorbed de- 
pending on our line-of-sight to the hotter or cooler side 
of this compact structure. Unlike in AGN, however, a 
dusty outflow associated with the puffed-up rim has not 
been observed. This difference is not unexpected given 
the much softer radiation field of a star compared to an 
AGN and the associated lower strength of UV/optical 
radiation pressure to launch such a wind. 

The question remains if objects like NGC 424 and 
Circinus (and possibly NGC 1068) are a special class of 
objects or if the dusty polar regions are more generic to 
the AGN population. As discussed above, there are signs 
of polar extensions in the mid-IR on scales of 10s of par- 
sees in other Seyfert galaxies as well, but we do not have 
the interferometry data to test if this polar component is 
energetically dominant in these sources. If it is generic, 
then it will require adjustments to current models of the 
dust emission of AGN. 

6. SUMMARY AND CONCLUSIONS 

We present VLTI/MIDI mid-IR interferometric obser- 
vations of the X-ray-obscured type 2 nucleus of NGC 424. 
The ww-coverage allowed us to constrain the position 
angle-dependent sizes of the mid-IR emission source in 
the 8 — 13 /im wavelength range. We used two differ- 
ent methods to overcome size biases by the inhomoge- 
neous distribution of projected baseline lengths on the 
uv-plane (i.e. beam shape effects). As complementary 
information, a high spatial resolution IR SED has been 
compiled. In addition to simple geometric fits, we also 
modeled photometric and interferometric data simulta- 
neously by a simple model that assumes dust emission 
in thermal equilibrium with the AGN. This led to the 
following main results: 

• The mid-IR emission source in NGC 424 is resolved 
and elongated. Its size depends on position-angle, 
wavelength, and baseline. For a baseline length of 
58m we find a size of (1.7±0.2) pcx (1.2±0.2) pc at 
9.1 //m, and(2.2±0.4) pcx(1.9±0.4) peat 12.5 /Ltm. 
As shown in Fig. [H the source size shows a system- 
atic increase with wavelength. A more baseline- 
independent size can be extracted usi ng the half- 
light radius R1/2 as defined in .Kishimoto et all 
(l2011bD . We obtain = (2.0 ± 0.2) pc x (1.5 ± 
0.3) pc for the 9—13 /im wavelength-averaged vis- 
ibilities of NGC 424. 

• The major axis of the mid-IR emission in NGC 424 



points toward PA —27° ± 2° based on the best-fit 
model of the complete data set. This direction is 
consistent with the PA derived from simple geomet- 
ric fits to individual wavelength data corrected for 
the uu-coverage bias (or beam shape) . We also de- 
tect a flatter radial brightness distribution in the 
direction of the major axis than along the minor 
axis of the mid-IR emission. Interestingly, the po- 
sition angle of the major axis is only 20° off the sys- 
tem axis as set by spectro-polarimetry. Therefore, 
we conclude that the mid-IR emission is extended 
in polar direction or symmetry axis of the torus. 

• The polar-elongated emission is responsible for 
>60% of the mid-IR total flux (i.e. single tele- 
scope flux) in NGC 424 on scales from about 1 pc 
to -100 pc. 

These results can be put into the broader context of 
AGN uniflcation and our current picture of the mid-IR 
emission from the dusty torus. We conclude: 

• The elongation of the nuclear mid-IR source on par- 
sec scales in polar direction seems inconsistent with 
typical state-of-the-art radiative transfer models of 
the torus, either smooth or cl umpy. While the 
hydro dynamic simulations by iSchartmann et all 
(|2009[ ) do show polar-elongated mid-IR emission, 
the model SEDs are incompatible with the ob- 
served IR SED of NGC 424. We conclude that it 
is very difficult to simultaneously explain both the 
shape and SED of NGC 424 within the framework 
of current torus models. 

• Based on the spectral (SED) and spatial (interfer- 
ometry) information, the source of the emission is 
most likely optically thin dust in the polar region of 
the AGN. While the total dust column and surface 
brightness are, therefore, low, the emission covers 
a relatively large area leading to the high contribu- 
tion of flux to the overall mid-IR emission. We note 
that similar mid-IR extensions in polar region have 
been observed on scales of several 10s of parsecs in 
a few other sources, including NGC 1068. 

• The characteristics of the mid-IR emission in 
NGC 424 show similarities to the nucleus of the 
Circinus galaxy (and probably NGC 1068 as well) 
where about ~85% of the mid-IR flux originates 
in the polar region and coincides with the western 
edge of the outflow cone., i.e. several degrees off 
from the nominal system axis. In NGC 424, the 
mid-IR major axis is off about 20° from the polar 
axis. This could be interpreted as an edge bright- 
ening effect in a cone filled with optically thin dust. 

• Our data suggests that the origin of the dust in 
the polar region may be a radiatively-driven wind 
from the inner part of the torus. The dust in the 
inner torus region is expected to be deficient of sil- 
icate grains, which is well matched by the mid-IR 
spectral features of the polar dust. Based on re- 
cent near- and mid-IR interferometry results of a 
sample of type 1 AGN, we propose that the torus 
has a small scale height at large distances from the 
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AGN and is puffed-up by radiation pressure in the 
inner region, leading to the required covering fac- 
tors in the unification scheme. In this picture, the 
near-IR emission (up to about 5 ^m would be dom- 
inated by the torus, while the dusty outflow would 
be the main contributor to the mid-IR. 



Our results for NGC 424 have significant implications 
for models of the IR emission of AGN. Owing to a lack 
of suitable data, we do not know if the mid-IR emis- 
sion sources in NGC 424 and a few other nearby Seyfert 
galaxies are special cases or more generic to AGN, but it 
seems as more cases of polar-extended mid-IR emission 
are revealed when objects are studied in detail by high- 
angular resolution tech niques. Moreove r , the r esults of 
lH5nig et all (|2010D and iKishimoto et all (|2011bD suggest 
that the structure of the dusty environment depends on 
luminosity. It is, therefore, important to follow-up our 
results with interferometry from near- to mid-IR wave- 
lengths and a good position-angle and baseline coverage 
of more sources. The upcoming VLTI/MATISSE instru- 
ment, operating from 3.5 to 13 ^m, should be able to 
test if the near-IR emission is, indeed, dominated by the 
torus while the mid-IR flux originates in the outflow. 
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APPENDIX 

NGC 424 VLTI/MIDI INTERFEROMETRY DATA (ONLINE-ONLY) 

In Fig. [TU] we present all reduced and calibrated correlated flux data (blue lines with error bars) and visibilities 
(gray lines) obtained from 2007 through 2010. The telescope configuration, projected baseline length, position angle, 
and observing dates are indicated. For more details on the data quality we refer to the maximum SNR values of 
each set indicated in Table [TJ As discussed in Sect. 14.2] the shortest wavelength data points are probably affected by 
decoherencing effects or the wing of the remaining 7.7 /im PAH feature. They have been excluded from the model 
fitting as well as the 2007 data points (see Sect. |4] for details). The best-fit model visibilities are overplotted as green- 
dashed lines in Fig. [TOl Note that some of the data sets show residuals of the atmospheric ozone feature at 9.6 ^m 
that remained after calibrating the data. 
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Figure 10. All 37 reduced and calibrated interfcrometric data sets obtained for NGC 424. The blue lines with error bars show the 
correlated fluxes. The visibilities are shown as gray lines with error bars. We also overplotted the binned visibilities that have been used 
for size extraction and model fitting (black-filled circles with error bars). The green-dashed line represents the hot ring -|- T/d-gradient 
model fit to all interfcrometric and photometric data as detailed in Sect. |4] 
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Figure 10. continued. 



